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Abstract 

The study of quark jets in e + e~ reactions at LEP has demonstrated that the 
hadronisation process is reproduced well by the Lund string model. However, 
our understanding of gluon fragmentation is less complete. In this study en- 
riched quark and gluon jet samples of different purities are selected in three-jet 
events from hadronic decays of the Z collected by the DELPHI experiment in 
the LEP runs during 1994 and 1995. The leading systems of the two kinds of 
jets are defined by requiring a rapidity gap and their sum of charges is studied. 
An excess of leading systems with total charge zero is found for gluon jets in 
all cases, when compared to Monte Carlo Simulations with JETSET (with and 
without Bose-Einstein correlations included) and ARIADNE. The correspond- 
ing leading systems of quark jets do not exhibit such an excess. The influence of 
the gap size and of the gluon purity on the effect is studied and a concentration 
of the excess of neutral leading systems at low invariant masses (< 2 GeV/c 2 ) 
is observed, indicating that gluon jets might have an additional hitherto unde- 
tected fragmentation mode via a two-gluon system. This could be an indication 
of a possible production of gluonic states as predicted by QCD. 
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1 Introduction 

The study of quark jets provides us with remarkable insights into the mechanism of 
hadronisation. It gives strong evidence for chain-like charge ordered particle production in 
excellent agreement with string Monte Carlo models like JETSET [1] . This is shown e.g. 
by several contributions [2-5] of the DELPHI experiment at LEP, where the compensation 
of quantum numbers, in particular that of charge, has been extensively studied. Much 
less is, however, known about the behaviour of gluon jets. On the theoretical side, besides 
the fragmentation via two strings as implemented in JETSET/PYTHIA and ARIADNE 
[6], the direct neutralisation of the color octet field by another gluon with the creation 
of a two-gluon system has been considered by Minkowski and Ochs [7,8] and also by 
Spiesberger and Zerwas [9]. Older references exist by Montvay [10] and Peterson and 
Walsh [11]. Additional references can be found in [7] where it is also emphasized that 
an experimental study of the gluon corner in three-jet events could contribute valuably 
to the question of the existence of glueballs, an early expectation of QCD [12]. No 
quantitative prediction however exists up to now. This has triggered an experimental 
investigation by the DELPHI collaboration on gluon fragmentation in a leading system 
defined by a rapidity gap [13,14]. The preliminary results revealed that electrically neutral 
systems of leading particles in gluon jets occur more often than predicted by JETSET, 
in agreement with the expectations of the above theoretical argumentations, while there 
was no disagreement observed in quark jets. This phenomenon, experimentally observed 
for the first time, has meanwhile also been seen by ALEPH and OPAL [15,16] 

The JETSET (ARIADNE) model of a qqg event stretches a string from the q to the 
g and on to the q. The string fragments for example by the creation of qq pairs, similar 
to what happens for quark fragmentation (Fig. la). Thus the JETSET (and ARIADNE) 
model regards gluon fragmentation as a double color triplet fragmentation (most clearly 
sketched in Fig.l of ref. [7]) and the leading system can obtain the charge ±1 or in 
the limiting configuration. The process proposed by Minkowski and Ochs, namely the 
octet neutralisation of the gluon field by another gluon has the signature of an uncharged 
leading system due to the requirement that the sum of charges (SQ) of the decay prod- 
ucts of a two-gluon system is zero (Fig. lb). In [7,8] it is also proposed to enhance the 
possible contribution of this process by selecting events where a leading particle system 
is separated from the rest of the low energy particles by a large rapidity gap, empty of 
hadrons. In this situation of a hard isolated gluon the octet field is expected not to have 
been distorted by multiple gluon emission and by related color neutralisation processes 
of small rapidity ranges [7]. The price to pay for such a selection is, however, a strong 
reduction of the number of events because of the Sudakov form factor [17]. A differ- 
ent mechanism - color reconnection [18] - can produce similar effects. Two experiments, 
however, agree that the present color reconnection models, as implemented in some ver- 
sions of Monte Carlo simulations, can not reproduce quantitatively the observed excess 
of SQ = systems [15,19]. 

The present study aims to consolidate the results of the preceding analyses [13-16] by 
studying the dependence of the excess of neutral leading systems in enriched gluon jets 
on the gap size and gluon content, by studying their mass spectra and by investigating 
if there are possible trivial origins for the observed effect. This is especially important, 
since a significant failure of the string model to describe gluon jets might generally reveal 
the presence of hitherto undetected processes. The size of the effect for a pure gluon jet 
is estimated. As a cross-check, the same investigation is done for quark jets. 
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Figure 1: Diagrams to illustrate the processes of color triplet fragmentation a) and color 
octet fragmentation b). The dashed lines represent the color triplet strings and the helixes 
represent the color octet strings. 

2 Data sample and 3-jet event selection 

The data sample used has been collected by the DELPHI experiment at the LEP 
collider at the peak of the Z resonance during 1994 and 1995. Three-jet events have 
been selected by using the appropriate cuts for track quality and for the hadronic event 
type [20] as well as applying a k t cluster algorithm (Durham) [21] to all observed charged 
and neutral particles with y cut = 0.015 1 . The jet energies were recalculated based on 
the direction of the jet momenta and the jets were sorted by decreasing energy, i.e. 
E3 < E 2 < Ei. Events with G 2 , 03 = 135° ±35° have been used, where 0j is the inter-jet 
angle opposite to jet i. All jets are required to lie in a plane, to consist of > 2 particles 
and the jets must be at least 30° away from the beam direction [22-24]. About 314000 
events meet all these conditions. 

Without any additional tag the jet with the highest energy E\ (jetl) is in most cases 
a quark jet and that with the smallest energy E3 (jet3) the gluon jet. The measured 
mean jet energies are: E x = 41 .4 GeV, E 2 = 32.2 GeV and E 3 = 17.7 GeV. 2 In the first 
data sample (samplel), where the gluon and quark jet identification is based on energy 
ordering only, events are required not to exhibit any b-signal (235080 events). Monte 
Carlo simulations show for the above mentioned conditions a quark jet contribution of 
about 90% for jetl 3 and a gluon jet contribution of about 70% for jet3. In a more detailed 
study of the gluon purity a second independent sample (sample2) is selected, where jetl 

1 This value has been obtained from a study optimizing simultaneously purity and statistics [22] 

2 Although the mean energies of jetl and jet3 differ by more than a factor 2, the maximum possible rapidities and mean 
multiplicities for charged particles differ much less (e.g. (nchj et :j) = 6.04, (nchj e ti) =7.36). 
3 A11 quark jet selections (jetl) shown in the figures for comparisons are defined by samplel. 
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and jet2, contrary to jet3, are required to exhibit a b-signal [24,25] (Section 4.3). This 
additional tag results in a gluon purity of jet3 of about 90% and consists of 31400 events. 
A third sample (sample3) is selected to enable purity unfolding for special cases. It is 
defined by the requirement that jet3 has a b-tag. For this jet3 sample consisting of 12200 
events the gluon content is very much diminished (about 26%). 

3 Monte Carlo models 

For comparisons a suitable number of Monte Carlo simulations using JETSET 7.3 [1] 
and ARIADNE [6] have been performed. In contrast to JETSET, ARIADNE incorpo- 
rates dipole radiation of gluons instead of the parton shower used by JETSET. Since 
Bose-Einstein correlations (BEC) are present in nature, like-charged particles will stick 
together in momentum space and local charge compensation is expected to be diminished. 
The implementation of Bose-Einstein correlations into the Monte Carlo simulation, how- 
ever, is highly problematic and the magnitude of the effect on charge compensation is 
unknown. Nevertheless, the possible effect of BEC has to be investigated and the possible 
uncertainties have to be considered. 

Three different Monte Carlo event samples have been created by using different gen- 
erators: 

Model (1): JETSET with BEC included (BE32 [26]) 
Model (2): JETSET without BEC 
Model (3): ARIADNE without BEC 

The number of events generated for each sample corresponds roughly to that of the 
data. 

The data are compared to these Monte Carlo event samples with full simulation of the 
DELPHI detector. The same reconstruction and analysis chain has been applied to the 
data and Monte Carlo (MC) samples. 

4 Analysis 

4.1 The sum of charges in the leading system with a rapidity 
gap (samplel) 

After the selection of 3-jet events and the determination of enriched quark and gluon 
jet samples, the leading system of a jet is defined by requiring that all charged particles 
assigned to the jet must have a rapidity y with respect to the jet axis of y > Ay. The 
quantity Ay represents a lower limit and defines a rapidity gap extending at least up 
to y = Ay. By this requirement, also jets are discarded, if they include some particles 
(fraction 10~ 3 ) with negative rapidity. The size of the demanded gap below this leading 
system is a compromise between the requirement of a gap as large as possible and the 
considerable loss of statistics at a larger gap. The requirement that the rapidity interval 
5y > Ay (with Ay =1.5) below the leading system be empty of charged particles reduces 
the number of jets appreciably. About 38000 enriched gluon jets and 39000 quark jets 
meet this condition. This reduction rate f± of enriched gluon jets is quite well reproduced 
by the three Monte Carlo event samples: /i(data) = 0.160, /i(MCl) = 0.169, /i(MC2) 
= 0.158, /i(MC3) = 0.157, with the mean value /i(MC-mean) = 0.161. In principle, 
there could be neutral hadrons in the gap. It has been verified that removing in addition 
topologies where observed neutrals are contained in the gap (mainly 7's from the decay 
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of 7r°'s), leads to results that are fully consistent with the ones presented here, but with 
about 15% larger statistical errors. 

The sum of charges (SQ) of the particles belonging to the leading system defined as 
above is shown in Fig. 2a for enriched gluon jets and in Fig.2b for quark jets and compared 
to ARIADNE. P{SQ, Ay) is generally defined as the fraction of a jet sample with a gap 
and a given value of SQ, 

and is an estimate for the probability of a jet with a gap to have a certain SQ. The 
SQ distribution of the leading system for the gluon jet (Fig. 2a) exhibits for SQ = a 
significant enhancement of the data over the Monte Carlo. This effect is predicted, if 
the process of color octet neutralisation is present [7,8]. On the other hand, there is no 
significant difference between the data and the Monte Carlo simulation in the case of 
quark jets (Fig.2b). 

The lower parts of Fig.2 show quantitatively the differences of the P(SQ, 1.5) between 
the data and the Monte Carlo simulation. This difference for the gluon jet (Fig. 2c) 
amounts to about 4 standard deviations (statistical errors only), for the quark jet (Fig.2d) 
this difference is compatible with zero. 

4.2 The dependence on the size of the rapidity gap 

Figs.3a,b show, for neutral leading systems (SQ = 0), the dependence of the relative 
deviation R(Ay) on the size of the lower limit (Ay) of the rapidity gaps considered 4 : 

= P( ( KA V) , P( (I ,A V)MC 
P(0,Ay) M c 

For all three types of Monte Carlo simulations (models (1), (2) and (3), see Section 3), 
R(Ay) (Ay > 0.5) is positive and increasing with Ay for jet3 (Fig.3a). This clearly shows 
that the surplus of neutral leading systems in the data, compared to the Monte Carlo 
simulations, increases with the gap size. This corroborates the arguments of Minkowski 
and Ochs [7,8]. ARIADNE without BEC (model 3) lies between the JETSET models. 
In the case of jetl (Fig.3b) all values are scattered around zero and no rise can be seen. 
When comparing JETSET with and without BEC included, one notices for jet3 and also 
for jetl a difference for all values of Ay. The effect of introducing BEC into the Monte 
Carlo models causes a shift of R(Ay), essentially independent of Ay, to higher values. 
The dependence on Ay is approximately the same for all three models. Arguing that only 
the rise and not implementation effects of BEC are of interest here and that the surplus 
of neutral systems is expected to be small and in a first approximation negligible in jets 
without a gap requirement, the following quantity is calculated: 

R' (Ay) = R(Ay) - R(0) (3) 

and shown in Fig. 3c, d. The residual spread between the models is considered as system- 
atic uncertainty. 

4.3 The dependence of R'(Ay) on the gluon purity 

In Fig. 4a the mean values of R' with models (1), (2) and (3), which are presented 
in Fig.3c for samplel, are drawn together with the statistical errors (symbols with error 

4 In this representation, the bins arc not independent: each point represents a subsamplc of the previous one. 
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Figure 2: Fraction of jets P(SQ, 1.5) as a function of the sum of charges SQ of the leading 
system for both (a) enriched gluon jets and (b) quark jets. Full circles represent the data, 
lines the Monte Carlo simulation, model (3). The difference (data-MC) is shown in (c) 
for enriched gluon jets and in (d) for quark jets. 
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Figure 3: a),b): Relative deviations R(Ay) of the number of neutral leading systems 
in gluon and quark jets for the three Monte Carlo models defined in Section 3 and for 
various sizes of Ay. R is defined in equ(2). c),d): R (Ay) as defined in equ(3). Because 
of the nature of the cut Ay, the bins are correlated. 
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Figure 4: R (Ay) using model averages. The gluon jet3 is tagged using energy ordering 
in a) and b-tagging is used for the selection in b). 

bars) and systematic uncertainty (shaded area). As a cross-check, an independent second 
sample (sample2) of gluon jets with a much higher purity is selected. The dependence 
of R (Ay) on Ay for this sample is given in Fig. 4b. Although the statistics are smaller 
in Fig.4b (3870 jet3 at Ay = 1.5, which is only about 1/10 of samplel), the effect is 
increased, which is expected if it is connected to the gluon jet only. At Ay = 1.5 R (Ay) 
is about 0.09 ± 0.02 (statistical) because of the higher purity. 

To estimate the amount of disagreement between data and Monte Carlo in pure gluon 
jets the gluon purity has to be estimated for jet3 at the gap size Ay for both data 
selections in Fig. 4. In principle, it can be directly obtained from the Monte Carlo. At the 
same scale, gluon and quark jets exhibit different rapidity distributions, i.e. gluon jets 
emit more particles per unit at small rapidity. Demanding a gap reduces therefore not 
only the number of jets, but also the gluon content in a mixed sample of gluon and quark 
jets. This is observed in the MC. An estimation of the gluon purity at gap Ay however 
depends on the correct modelling of the rapidity distribution of pure gluon jets. Therefore 
another method has been used in addition. It uses the measured reduction rates of the 
number of jets U(Ay) in sample-i by demanding a gap (see Section 4.1) and from the MC 
only the composition at gap zero. Let us define N^Ay) = / 1 (Ay)iV 1 (0) in samplel, and 
N 2 (Ay) = f 2 (Ay)N 2 (0) in sample2, where Ni(Ay) (N 2 (Ay)) is the number of jets counted 
at Ay in samplel (sample2). Since these samples are an admixture of pure q(=light quark) 
jets, g(=gluon) jets and b(=b-quark) jets, the corresponding fi(Ay),f 2 (Ay) (and also 
f 3 (Ay) for sample3) are also an admixture of the reduction rates f q (Ay), f g (Ay), fb(Ay) 
of the pure light-quark, pure gluon and pure b-quark subsamples, e.g. 

ft(Ay) = a lq f q (Ay) + a lg f g (Ay) + a lb f b (Ay) (4) 
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with two analogous equations for f 2 (Ay) and / 3 (Ay). The resulting system of three linear 
equations can be written in short: 

F = AF pure (5) 

with the solution 

F pure = A~ l F. (6) 

The vector F(f 1 (Ay), f 2 (Ay), / 3 (Ay)) is measured, 

and the vector F pure (f q (Ay), f g (Ay), f b (Ay)) is the solution. 

The matrix A represents the q,g,b compositions for the three selections at gap=0, es- 
timated from Monte Carlo (e.g. a\ g is the gluon purity of jet3 in sampler, a 2g that of 
sample2, and a^ g that of sample3 and so on). With the solution of equ.(6), the numbers 
of true gluon jets at Ay can be determined in sample 1 and sample2: 

Nt on (Ay) = f g (Ay)a lg N l (0), (7) 

Nt on {Ay) = f g (Ay)a 2g N 2 (0). (8) 
The fraction of gluon jets c^ y at Ay is given by: 

cfy(samplel) = Nf luon (Ay)/N 1 (Ay) = a lg f g {Ay)/h{Ay), (9) 

cfy(sample2) = N? uon (Ay) / N 2 (Ay) = a 2g f g (Ay)/f 2 (Ay). (10) 

Applied to the two data sets in Fig.4a,b the following numbers for the gluon content 
have been obtained: 

1. The sample in Fig. 4a (sample 1): c Q g = 0.65, c*' 5 = 0.46 (from equ.(9) ), and c l g b = 
0.45 (directly from the Monte Carlo at Ay = 1.5). 

2. The sample in Fig.4b (sample 2): c° g = 0.88, c^ 5 = 0.80 (from equ.(10) ), and 
cj- 5 = 0.82 (directly from the Monte Carlo at Ay = 1.5). 

The statistical errors on these numbers are below 1%, systematic errors can be obtained 
by comparing the estimates with different Monte Carlo models (1), (2) and (3). They 
are < 2.6%. The purity estimates obtained above allow the determination of the excess 
of neutral systems R g in pure gluon jets by dividing by c^' 5 . The following values of R g 
have been obtained for the two samples defined above: 

Samplel : # fl (l-5) = 0.100 ± 0.012 (stat) ± 0.025 (syst) (11) 

Sample2 : i^(1.5) = 0.107 ± 0.022 (stat) ± 0.028 (syst) (12) 

The samples are statistically independent. Adding statistical and systematic errors 
quadratically, a significance of about 3.6<r is obtained in sample 1 and of almost 3<r 
in sample 2. Combining finally both samples results in: 

Combined : ^(1-5) = 0.102 ± 0.011 (stat) ± 0.026 (syst) (13) 

This number can be used to make a first estimate of R g (0), the amount of the excess of 
neutral systems in pure gluon jets without any gap selection. Taking into account the 
estimated value of / ff (1.5) = 0.112±0.003 from equ.(6) which tells that 11% of the pure 
gluon jets meet the gap condition Ay = 1.5, one obtains: 



R g (0) ~ R' g (Ay = 1.5)f g (Ay = 1.5) = 0.011 



(14) 
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Extending this analysis to samples which allow also for smaller gap sizes 1.5 > Ay >1.0 
leads to R g (0) values between 0.01 and 0.02. Extending this furthermore to all gluon jets 
by taking into account that P 9 (0,0) ~ 0.26 (from sample2, not shown here) leads to the 
conclusion that the amount of a possible octet neutralisation of the gluon field is of the 
order of 0.5%. 

4.4 Discussion of the systematic uncertainties 

The following sources of systematic errors have been considered: 

(a) Quality of event reconstruction. Bad reconstructions and losses of particles (mainly 
neutrals) in the detector and wrong assignments to the jets can lead to differences of 
several GeV between the jet energy calculated from the angles between jets (E ca [ c ) 
[24] and the sum of energies of all particles assigned to the jet (E sum ). Improving the 
quality by cutting away about 1/3 of the jets with the largest difference E ca i c — E sum 
does not significantly change the signals at SQ = both in gluon and quark jets. 

(b) The dependence of the effect on the polar angle of the jet with respect to the collision 
axis has been investigated: the effect is stable. 

(c) The influence of track finding efficiency in the detector. In order to investigate 
the influence of track finding efficiency the effect of a reduction of the efficiency by 
1% has been simulated. No significant change in the signals at SQ = has been 
observed. Since R is a ratio with respect to the Monte Carlo simulation (including 
the detector effects) and the deviation which has been observed between data and 
MC below 10%, it can be expected that efficiency effects cancel to a large extent. 

(d) To investigate whether the good agreement between data and Monte Carlo in quark 
jets is only due to the larger particle momenta, in a test-run only particles with 
momenta less than 30 GeV/c have been accepted in jetl. The agreement with the 
Monte Carlo remains. 

(e) The estimations leading to (11) and (12) assume that quark jets, also at the lower 
energies of jet3, do not exhibit any excess of neutral leading systems. This is further 
tested by measuring the excess in sample3 which exhibits at Ay = 1.5 an admixture 
of only 20% gluon jets. As expected, the signal is reduced, and is even negative with 
large error: #3(1.5) = -0.02 ± 0.02. Adopting the same procedure as in Section 
4.3 by using matrix inversion with the measured values R^l.5), i= 1,2,3 for the 3 
selected samples as input, the resulting excess for pure quark and gluon jets could 
be estimated: R' q (1.5) = 0.00 ± 0.02, R' g (1.5) = 0.11 ± 0.03 and R' b {1.5) = -0.06 ± 
0.04. These results do not show any evidence that quark jets exhibit an excess of 
neutral leading systems for the lower jet3 energies. 

(f) At the generator level of JETSET and for pure gluon jets the effect of changing 
parameters within limits has been studied. For example, different DELPHI tunings 
have been used, the DELPHI tuning [27] has been replaced by that of OPAL [28] 
and by the JETSET default 5 , and the popcorn parameter has been varied. Some 
changes of P(SQ, Ay = 1.5) at SQ = are revealed in gluon jets and to a lesser 
extent in pure quark jets. At a gap size of Ay = 1.5 a maximum variation of R(Ay) 
of 0.027 is observed. 

The systematic error from (f) amounting to 18% is taken into account. This is a 
conservative estimate with a factor 0.68 of the maximum variation, corresponding to la 
of a Gauss distribution. The contributions from (a)-(e) are negligible. The systematic 
errors for samples 1 and 2 are estimated as follows: 

5 All these studies have been done with BE correlations included. 
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1) Sample 1 

a) from the spread in Fig. 4a at Ay = 1.5: 
AR' = 0.017 

b) uncertainty in purity: 0.0026 (see Section 4.3) 

c) uncertainty from (f): 0.018 

By quadratically adding all 3 contributions a systematic error of 0.025 is obtained. 

2) Sample 2 

a) from the spread in Fig. 4b at Ay = 1.5: 
AR' = 0.021 

b) uncertainty in purity: 0.0028 (see Section 4.3) 

c) uncertainty from (f): 0.019 

By quadratically adding all 3 contributions a systematic error of 0.028 is obtained. 

5 Mass spectra 

Color octet neutralisation of the gluon field could produce a resonance spectrum which 
differs from that of color triplet fragmentation [29] implemented in JETSET. In order 
to investigate in which region of the mass spectrum the observed excess of the leading 
neutral systems is located, the invariant mass (M) distributions 

^m^M (18) 

of the leading systems with total charge zero at Ay = 1.5 have been calculated and 
compared with the mean values of models 1-3 for the two cases: 

a. M is computed using only charged particles with a momentum p > 0.2 GeV/c and 
assuming pion mass. This distribution is shown in Fig.5a for gluon enriched jets 
and in Fig. 5c for quark jets. Two broad bumps can be observed in Figs. 5a, c which 
are the result of a superposition of a rapidly decreasing two-particle spectrum which 
dominates for M < 1 GeV/c 2 , with an increasing spectrum consisting of 4 and more 
particles. The latter dominates and peaks at ~ 1.5 GeV/c 2 . The region below 0.8 
GeV/c 2 consists only of two particles. One peak around M ~ 0.8 GeV/c 2 can be 
attributed to the p resonance, another at M < 0.5 GeV/c 2 to a reflection of i], rj' and 
uj. The latter statement is corroborated by the fact that in events with no neutrals, 
the peak at M < 0.5 GeV/c 2 vanishes. Only in Fig.5a a third peak is indicated by 
the data points just below 1 GeV/c 2 , in the region of the /o(980) resonance. In this 
region 0.9 < M < 1 GeV/c 2 the two particle contribution amounts to about 70%. 
Figs. 5b, d show the difference between the data and the Monte Carlo event sample. 
For gluon enriched jets the distribution in Fig.5b exhibits possible evidence for a 
mass enhancement 6 in the region 0.9 < M < 0.95 GeV/c 2 , but no signal is seen in 
this region for quark jets in Fig.5d. This peak is the remaining part of the original 
peak in Fig. 5a after the subtraction of a small and narrow but significant signal in 
the MC at 0.95-1.0 GeV/c 2 . Without emphas izing too much this remaining narrow 
peak in Fig.5b, it has to be noted that it survived a quality cut (by accepting only 
jets with a polar angle > 50°) well above 3a, whereas all other deviations from zero 
in isolated bins in Figs. 5b, d were decreased. Whether this narrow signal in gluon 

6 The experimental mass resolution is below 10 MeV/c 2 , the binwidth in Fig. 5 is 50 MeV/c 2 . 
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Figure 5: a),c): Invariant mass distribution P(M) of the leading system (SQ = 0) 
(considering only charged particles) a)for gluon-enriched jets (sampler), c) for quark jets. 
P is defined in equ.(15). The dots with error bars are the data, the histograms are the 
mean values of the three Monte Carlos. b),d): P data (M) - P MC (M). 



enriched jets can be attributed to / (980) production remains an open interesting 
question for future studies with increased statistics. In the mass range below 2 
GeV/c 2 an overall excess is observed in Fig. 5b. 

b. M is computed for all charged particles (with p > 0.2 GeV/ c and pion mass assumed) 
and neutrals (photons with E > 0.5 GeV). Figs.6a,d show P(M) for gluon enriched 
jets and Fig.6g for quark jets. Both gluon enriched samples exhibit an excess of 
data compared to the MC for low invariant masses (below 2 GeV/c 2 ), which is 
emphasized in Figs.6b,e where the mass distributions of data and MC have been 
subtracted. Separating the mass ranges below and above 2 GeV/c 2 in Fig. 6c for 
sampler, an excess of about 5<r(stat) is observed in the low mass range. Comparing 
Fig. 6c with Fig.6f, one observes that the excess is increasing according to increasing 
gluon purity, namely about a factor 2 between samplel and sample2. For quark jets, 
the corresponding distributions (Fig.6g,h,i) do not exhibit any significant difference 
between data and MC. 

An important remark concerns all mass spectra in Fig. 6. As stated in section 3, all 
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the comparisons to Monte Carlo event samples are done with full detector simulation. 
Consequently, due to the loss of neutral particles, all spectra and in particular the 
excess spectra in Figs. 6b, 6e are shifted by about 0.3 to 0.5 GeV/c 2 to lower mass 
values. It has been verified, however, by a special Monte Carlo study using the 
detector response matrix (not shown here) that, after the correction of all spectra 
for this shift, the excess is still clearly concentrated in the low mass region. On 
the other hand, the spectra of leading systems consisting only of charged particles 
(Fig. 5) are not affected by shifts. 

The observations in Fig. 5 confirm the first preliminary results presented in 2001 [13] 
for leading gluonic systems by considering charged particles only. In 2002 the OPAL 
collaboration reported [16] a 2a excess in the mass distribution of neutral leading systems, 
consisting of charged and neutral particles, between 1 and 2.5 GeV/c 2 in gluon jets. 

The observation that the excess of neutral leading systems in gluon jets is limited to 
the low mass region (Figs. 5 and 6) supports arguments in favour of gluonic states. The 
existence of glueballs, i.e. bound states of two or more gluons, is a prediction of QCD 
[12]. The experimental results and their interpretations, however, are still controversial 
[30]. Theoretically there is general agreement that the lightest glueball should be in 
the scalar channel with J PC = ++ . Quantitative results are derived from the QCD 
lattice calculations [31] which predict the lightest glueball to be around 1600 MeV/c 2 , 
or from QCD sum rules [32], which predict also a possible gluonic state near 1 GeV/c 2 . 
Alternatively, it could also be a very broad object [33,34]. The gluonic state could mix 
with ordinary 0++ states, like the / (980), /o(1370), / (1500) or / (1710). As an example 
there are scenarios, where the largest gluonic component is included in the /o(1500) [35], 
or alternatively in the /o(1710) [36]. Recent discussions with various references can be 
found in [30,37]. 
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Figure 6: a),d): Invariant mass distributions P(M) of the leading system (SQ = 0) (with 
charged and neutral particles) for gluon-enriched jets, samples 1 and 2 respectively, g): 
same for quark jets (jetl, sampler). The dots with error bars are the data, the histograms 
are the mean values of the three Monte Carlos. The quantity P(M) is defined in equ.(15). 
Second row b),e),h): P data (M) - P MC (M); Last row c),f),i): P data (M) - P MC (M) for 2 
bins: (0.25-2 GeV/c 2 ) and (2-3.75 GeV/c 2 ). 
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6 Summary 

In the present study the leading systems defined by a rapidity gap have been inves- 
tigated for gluon and quark jets. The statistics of 1994 and 1995 at ^/s = 91.2 GeV 
obtained by the DELPHI collaboration is used to select 3-jet events and to single out 
quark jets (purity ~ 90%) and gluon enriched jets (purity ~ 70%) by energy ordering 
(samplel). For the (enriched) gluon jets a higher rate of neutral leading systems than pre- 
dicted by the Lund string model JETSET/ARIADNE (with and without Bose Einstein 
correlations) is observed but no such enhancement is seen for the quark jets. Various 
checks have been performed which suggest that this effect is not a spurious one. An 
increase of the effect with increasing gluon purity, obtained by a tagging procedure in a 
second sample (sample2), is observed corroborating that it is indeed connected with the 
gluon jets. 

The excess of neutral leading systems in pure gluon jets at a gap size Ay = 1.5 has 
been measured to be about 10%, with a significance of 3.6a. It is expected to be of the 
order of 0.5% in pure gluon jets without any charge or gap selection. 

The mass spectra of the neutral leading systems of gluon jets, both with and without 
including neutral particles have been studied. Mass spectra which include charged and 
neutral particles, show clearly that the excess mentioned above is concentrated at low 
invariant masses (< 2 GeV/c 2 ). The significance is enhanced there and amounts to about 
5a (statistical) in samplel and the excess is increased roughly proportionally to the gluon 
purity in sample2. 

The corresponding mass spectra of leading systems in quark jets do not exhibit any 
excess in the low mass regions. 

The observed excess of neutral systems in gluon jets and its increase with the gap size 
and with the gluon purity is in agreement with expectations, if the hitherto unobserved 
but predicted process of octet neutralisation of the gluon field takes place in nature. Al- 
though color reconnection could in principle alternatively explain the excess, the specific 
mass concentration at low mass seems to favor the first case and could be a signal of 
gluonic states predicted by QCD. 
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